How does degree heterogeneity affect nucleation of Ising model on complex networks? 
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We investigate the nucleation of Ising model on complex networks and focus on the role played by 
the heterogeneity of degree distribution on nucleation rate. Using Monte Carlo simulation combined 
with forward flux sampling, we find that for a weak external field the nucleation rate decreases 
monotonically as degree heterogeneity increases. Interestingly, for a relatively strong external field 
the nucleation rate exhibits a nonmonotonic dependence on degree heterogeneity, in which there 
exists a maximal nucleation rate at an intermediate level of degree heterogeneity. Furthermore, we 
develop a heterogeneous mean-field theory for evaluating the free-energy barrier of nucleation. The 
theoretical estimations are qualitatively consistent with the simulation results. Our study suggests 
that degree heterogeneity plays a nontrivial role in the dynamics of phase transition in networked 
Ising systems. 

PACS numbers: 89.75.Hc, 64.60.Q-, 05.50.+q 



I. INTRODUCTION 



Since many social, biological, and physical systems can 
be properly described by complex networks, dynamics on 
complex networks have received considerable attention in 
recent decades [H-Q- In particular, phase transitions on 
complex networks have been a subject of intense research 
in the field of statistical physics and many other disci- 
plines Q . Owing to the heterogeneity in degree distribu- 
tion, phase transitions on complex networks are drasti- 
cally different from those on regular lattices in Euclidean 
space. For instance, degree heterogeneity can lead to a 
vanishing percolation threshold @ , the whole infection of 
disease with any small spreading rate TjL the Ising model 
to be ordered at all temperatures [§ 4l0|, the transition 
from order to disorder in voter models [ll| , synchroniza- 
tion to be suppressed [H, EH and different path towards 
synchronization in oscillator network [3], just to list a 
few. However, there is much less attention paid to the 
dynamics of phase transition itself on complex networks, 
such as nucleation process in a first-order phase transi- 
tion. 

Nucleation is a fluctuation-driven process that initiates 
the decay of a metastable state into a more stable one 
[l5j . Many important phenomena in nature, like crys- 
tallization 16], glass formation [17] . and protein folding 
[HI are closely related to nucleation process. In the con- 
text of complex networks, the study of nucleation process 
is not only of theoretical importance for understanding 
how a first-order phase transition happens in networked 
systems, but also may have potential implications in real 
situations, such as the transitions between different dy- 
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namical attractors in neural networks and the ge- 
netic switch between high and low-expression states in 
gen e regulatory networks [2(J [HJ , and opinion revolution 
j22j as well as language replacement [23|, [24j in social 
networks. 

Recently, we have made the first step for studying 
nucleation process of Ising model on complex networks, 
where we have identified that nucleation pathways using 
rare-event sampling technique, such as nucleating from 
nodes with smaller degree on heterogeneous networks [25[ 
and multi-step nucleation process on modular networks 
[26| . In addition, we found that the size-effect of the nu- 
cleation rate on mean-field-type networks [1^] and non- 
monotonic dependence of the nucleation rate on modu- 
larity of networks (26|. As mentioned above, degree het- 
erogeneity has a significant effect on dynamics on com- 
plex networks. Therefore, a natural question arises: how 
degree heterogeneity affects nucleation of Ising model on 
complex networks? To answer this question, in this pa- 
per, we study the dynamics of nucleation on various net- 
work models whose heterogeneity of degree distribution 
can be continuously changed by adjusting a single pa- 
rameter. We use Monte Carlo (MC) simulation combined 
with forward flux sampling (FFS) to compute the nucle- 
ation rate and consider the effect of degree heterogeneity 
on the rate. Since the critical temperature of Ising model 
on uncorrelated random networks increases with the het- 
erogeneity of degree distribution 

(Emu, one may come 

to an intuitive conclusion: if both the temperature and 
external field are fixed, the nucleation rate will decrease 
monotonically as degree heterogeneity increases. Here, 
we show that such an intuition is not the case: the nucle- 
ation rate can change monotonically or nonmonotonically 
with degree heterogeneity depending on the level of driv- 
ing force, i.e., the value of external field. For a weak ex- 
ternal field, the nucleation rate decreases monotonically 
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with degree heterogeneity, whereas for a relatively strong 
external field there exists a maximal nucleation rate cor- 
responding to a moderate level of degree heterogeneity. 
Furthermore, we present a heterogeneous mean-field the- 
ory for calculating free-energy barrier of nucleation. The 
theoretical results qualitatively agree with the simulation 



II. MODEL AND SIMULATION 
DESCRIPTIONS 

The Ising model in a network comprised of N nodes is 
described by the Hamiltonian 



(1) 



t<3 



where spin variable Sj at node i takes either +1 (up) or 
— 1 (down). J(> 0) is the coupling constant and h is the 
external field imposed on each node. The elements of the 
adjacency matrix of the network take = 1 if nodes i 
and j are connected and a,j = otherwise. 

The simulation is performed by standard Metropolis 
spin-flip dynamics, in which we attempt to flip each spin 
once, on average, during each MC cycle. In each attempt, 
a randomly chosen spin is flipped with the probability 
min(l, e _/3AE ), where (3 = l/ffegJ 1 ) with the Boltzmann 
constant ks and the temperature T, and AE is the en- 
ergy change due to the flipping process. We set h > 
and T < T c , where T c is the critical temperature. The 
initial configuration is prepared with a metastable state 
in which Sj = — 1 for most of the spins. The system will 
stay in that state for a significantly long time before un- 
dergoing a nucleating transition to the thermodynamic 
stable state with most spins pointing up. 

Since nucleation is an activated process that occurs ex- 
tremely slow, brute-force simulation is prohibitively ex- 
pensive. To overcome this difficulty, we will use a recently 
developed simulation method, FFS [27j ■ This method al- 
lows us to calculate nucleation rate and determine the 
properties of ensemble toward nucleation pathways. The 
simulation results below are obtained by averaging over 
at least 5 independent FFS samplings and 10 different 
network realizations. 



III. RESULTS 

To study the effect of degree heterogeneity on nu- 
cleation, we first adopt a network model proposed in 
Ref.(28[. The network model allows us to construct net- 
works with the same mean degree, interpolating from 
Erdo-Rehyi (ER) graphs to Barabasi-Albert (BA) SF 
networks by tuning a single parameter 5 ebb A- For 
Serb A = one gets ER graphs with a Poissonian de- 
gree distribution whereas for Serb A = 1 the resulting 
networks are SF with P(k) ~ k~ 3 . Increasing Serb A 
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FIG. 1: The logarithm of the nucleation rate lni? as a func- 
tion of the strength of degree heterogeneity Serba for h — 0.5 
(a), h = 0.8 (b), and h — 1.0 (c). Other parameters are 
N = 1000, the mean degree (k) = 6, and T = 2.5. 



from to 1, the degree heterogeneity of the network in- 
creases. FigQ] shows that the logarithm of the nucle- 
ation rate lni? as a function of S for three different ex- 
ternal fields: h = 0.5, 0.8, and 1.0. For h = 0.5, lni? 
decreases monotonically with Serba, implying that de- 
gree heterogeneity is unfavorable for the occurrence of 
nucleation events. Interestingly, for h = 0.8 lni? is no 
longer monotonically dependent on Serba- as degree 
heterogeneity increases, lni? first increases slowly until 
Serba = 0.5 and then decreases rapidly. Further increas- 
ing h to h = 1.0, lni? clearly exhibits a nonmonotonic 
change with Serba- That is, there exists a maximal nu- 
cleation rate that occurs at a moderate strength of degree 
heterogeneity. 

To understand the above simulation results, we shall 
give a heterogenous mean-field theory on complex net- 
works for evaluating the nucleation barrier. First, we 
define rrik as the average magnetization of a node with 
degree k, i.e., m k = N^ 1 ^2 i ^ k _ k s i , where N k is the 
number of nodes with degree k. Furthermore, for a net- 
work without degree correlation, the probability that 
a randomly chosen nearest neighbor node has degree 
k is kP(k)/(k), where P(k) — Nk/N is degree dis- 
tribution and (k) = J^k kP(k) is the mean degree. 
Thus, the interaction energy between a node with de- 
gree k and its neighboring nodes can be expressed as 
—Jkrrik J2k' k' P(k')nik' / (k). The total energy of the net- 
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FIG. 2: Theoretical results of — f5AF as a function of Serba 
for h — 0.5 (a), h — 0.8 (b), and h = 1.0 (c). Other parame- 
ters are the same as those in Fig[T] 



work can be written as 



k'P(k')m k > 



(k) 



m k 



--NJ(k) to' 2 ~Nhm, 



where 



kP(k)m k 

W) 



(2) 



(3) 



is the average magnetization of a randomly chosen near- 
est neighbor node, and to = J^ fe P(k)m k is the aver- 
age magnetization of a randomly chosen node. Note 
that ml differs from to in general. Special cases for 
which to' = m are provided by fc-independent quantities 
TOfe = to. In particular, for the all-spin down configura- 
tion with TOfe = — 1 for all k and for the all-spin- up con- 
figuration with mfc = 1 for all k, one has to' = to = — 1 



and 



to = 1, respectively. 



Defining S k as the entropy of a node with degree k, 
the total entropy of the network is 



S = ^,N k S k = Nj2P{k)S k , 



(4) 



with 

S k = —Ub 



1 + rn k ^ / 1 + TO fc \ 1 - TO fc , n / 1 - TO fe 



FIG. 3: Simulation (left panels) and theoretical (right pan- 
els) results on networks with uniform degree distribution. The 
external fields from top to bottom are h = 1.0, 2.0, and 3.0, 
respectively. Other parameters are N = 1000, the mean de- 
gree (k) = 10, and T = 3. 



Combining Eqj2] and Eq|4j we can get the expression of 
free energy, i.e., F = E — TS. 

At the minimum and maximum points of free energy, 
we have dFjdm k = 0, which yields the mean-field equa- 
tion of m k [a, II , 



m k = tanh [j3h + [3 J km'] 
Substituting Eq(S]into EqJ31 we get 

fcP(fc) 



(k) 



■ tanh [(3h + ft J km'] 



(6) 



(7) 



(5) 



EqJT] is a self-consistent equation of to' that can be nu- 
merically solved. In the present settings, Eqj7]has three 
solutions: m'_, m' , and m',, where to' ± are stable solu- 
tions and to.q is unstable one. Inserting the three solu- 
tions of to' into the rsh of EqjBl we can obtain m k , and 
then get E a , S a and F a (a = — , 0, +) according to Eqj2] 
and EqHl Since h > 0, we have Fq > F- > F+, which 
gives the free-energy barrier from metastable to stable 
state AF = Fo — F- and thus estimate the nucleation 
rate R ~ exp(-/3AF). 

Theoretical results of —/3AF as a function of 5 are 
shown in Fig. 2, where the parameters are the same as 
those in Fig[TJ It is clear that the theoretical results are 
qualitatively consistent with the simulation ones. 

In order to check the generality of the above results, 
we shall calculate nucleation rate on some other net- 
work models by both numerical simulations and theory. 
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FIG. 4: Simulation (left panels) and theoretical (right pan- 
els) results on networks with Gaussian degree distribution. 
The external fields from top to bottom are h = 1.0, 2.0, and 
3.0, respectively. Other parameters are TV = 1000, the mean 
degree (k) = 10, and T = 3. 



Firstly, we construct a network with uniform degree dis- 
tribution in which node degree is randomly selected in 
the range [(k) - 5 mi , (k) + 5 uni ], where 5 uni is an integer 
between and (k) — 1 that controls the strength of degree 
heterogeneity. The network is generated according to the 
Molloy-Reed algorithm [2!| . This construction eliminates 
the degree correlations between neighboring nodes. Fig|3] 
shows the simulation and theoretical results, from which 
the similar phenomena are also present: for weak exter- 
nal field the nucleation rate decreases monotonically with 
degree heterogeneity, while for strong external field the 



nucleation rate varies nonmonotonically with degree het- 
erogeneity. Moreover, we construct a network with Gaus- 
sian degree distribution with fixed mean degree (k) and 
variance 6 gau . As shown in FigUl both the simulation 
and theoretical results display the similar phenomena. 



IV. SUMMARY 

In summary, using Ising model on complex networks 
we have shown how degree heterogeneity affects the rate 
of nucleation. The main results of the present paper 
are that for a weak external field the nucleation rate de- 
creases monotonically as degree heterogeneity increases, 
whereas for a relatively strong external field the nucle- 
ation rate first increases and then decreases with the 
increment of degree heterogeneity. Therefore, the nucle- 
ation rate can change monotonically or nonmonotonically 
with degree heterogeneity depending on the value of the 
external field. The results are robust to different network 
models, thereby verifying the generality of the results. 
Moreover, we have developed the so-called heterogeneous 
mean-field theory for calculating the free-energy barrier 
to nucleate and thus estimating the nucleation rate. The 
theory is effective in qualitatively predicting the simula- 
tion results. Our findings indicate that degree hetero- 
geneity plays a nontrivial role in the nucleation events of 
Ising model on complex networks. 
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